An ab-initio evaluation of the local effective interactions in the NaxCo02 family 
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We used quantum chemical ab initio methods to determine the effective parameters of Hubbard 
and t — J models for the NaxCo02 compounds (x=0 and 0.5). As for the superconducting compound 
we found the aig cobalt orbitals above the e'g ones by a few hundreds of meV due ti the e'g-eg 
hybridization of the cobalt 3d orbitals. The correlation strength was found to increase with the 
sodium content x while the in-plane AFM coupling decreases. The less correlated system was found 
to be the pure C0O2, however it is still strongly correlated and very close =to the Mott transition. 
Indeed we found U/t ^ 15, which is the critical value for the Mott transition in a triangular lattice. 
Finally, one finds the magnetic exchanges in the C0O2 layers, strongly dependant of the weak local 
structural distortions. 



I. INTRODUCTION 

Since the discovery of the large thermoelectric power^ 
and then of superconductivity^ in the NaxCo02 layered 
cobaltate family, these systems have attracted a lot of 
attention. Indeed, they present a very rich phases dia- 
gram, as a function of the sodium content x, and the pos- 
sible water intercalation. For small a: (0 < a; < 0.5) the 
magnetic susceptibility does not depend on temperature. 
These systems are thus supposed to be paramagnetic 
metals, hovi^ever some authors reinterpreted the data (us- 
ing NMR measurements^) and conclude that these sys- 
tems should be antiferromagnetic metals. For larger x 
(0.5 < a; < 0.75) the compounds are Curie- Weiss met- 
als. Over X = 0.75, the system has be seen as a weakly 
ferromagnetic mctal^, however neutron diffraction mea- 
surements^ rather sees A-typc antifcrromagnctism, that 
is ferromagnetic C0O2 layers, antiferromagnetically cou- 
pled. For special values of the sodium content, x = 1/4S., 
X = 1/2^ and x = 1—, the system is insulating. A band 
insulator for a; = 1 and a charge/spin ordered state for 
X — 0.5. On the contrary, for a; = the system remains 
conducting^^— , exhibiting a Fermi liquid behavior, while 
its half-filling character associated with the large electron 
correlation of the cobalt 3d let expect a Mott insulator. 

The NaxCo02 compounds are composed of C0O2 lay- 
ers in the (a, b) plane, between which the sodium ions 
are located. The C0O2 layers are formed by CoOe edge- 
sharing octahedra forming a two-dimensional triangular 
lattice of cobalt ions. The CoOg octahedra present a 
trigonal distortion along the crystallographic c axis cor- 
responding to a compression of the oxygen planes toward 
the cobalt one. Additional distortions take place, accord- 
ing to the sodium ordering. 

In the pure C0O2 compound, although Tarrascon et 
suggested a complex crystallographic structure, fur- 
ther authors refined the structure with a single cobalt 
crystallographic position -'9. This system exhibits a ge- 
ometry quite characteristic compared to the rest of the 
NaxCo02 family. First, the inter-layer distance is very 
short. Indeed, the inter-layer distance is 4.48 A , 4.31 A 
or 4.25A according to the authors, while for a; 7^ it 



ranges from 5.2 A to 5.7 A and increases with decreas- 
ing X. Second, the C0O2 layer is very compressed in the 
c direction. Indeed, its thickness is only 1.71 A while 
it ranges between 1.92 A and 1.98 A for the NaxCo02 
family. In fact, in this system, the C0O2 layer is even 
more compressed than in the superconducting, hydrated 
compound, for which its thickness is 1.76 A. It results in 
a very large trigonal distortion of the CoOe octahedra, 
with 62.3 degrees between the c axis and the Co-0 di- 
rections. Using the results of reference [l^, this strong 
trigonal distortion lets us expect a strong aig-e'g split- 
ting. Third, the Co-0 distances are shorter than in all 
the other compounds of the family, with only 1.847 A. 
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FIG. 1: Schematic structure of the Nao.5Co02 and NaoCo02 
compounds. The black lines connect the Na"*" ions with their 
nearest cobalt neighbors. 



In the Nao.5Co02 compound (P2 structure, subject of 
the present work), the sodiums ordering is commensu- 
rate and presents two types of alternated stripes along 
the b directioni^. In the first stripes, the sodium atoms 
are located on-top of (below) the cobalts atoms, and in 
the second type of stripes the sodium ions are located 
on-top of (below) the center of the cobalt triangles (see 
figure[T]). It results two types of cobalt atoms organized in 
stripes along the b direction, namely the Co(l) atoms (lo- 
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cated on top of a sodium ion) and the Co(2) ones. Com- 
pared to the other systems of the family, the Nao.5Co02 
compound does not stick on^. Indeed, its C0O2 layers 
present a thickness of 1.97 A . The average Co-0 dis- 
tances are 1.887A for Co(l) and 1.900A for Co(2). The 
average trigonal distortions are associated with angles of 
58.9 degrees (between the c axis and the Co(l)-0 di- 
rections), and 59.4 degrees (Co(2)-0). These angles are 
weaker than for the superconducting system or the pure 
C0O2 one, but still quite larger than the 54.7 degrees 
of the regular octahedron. One can thus expect a split- 
ting of the t2g cobalt 3d orbitals in two quasi-degenerate 
Cg-type orbitals of low energy and a aig-type one at the 
Fermi level. 



Out of all the NaxCo02 systems, the x — 0.5 com- 
pound presents a very rich but not well understood 
phases diagram. Indeed, the Na"*" ions order at very 
large temperature, namely slightly above 300 K^^, in- 
ducing a small charge order in the cobalt layer. Despite 
this small charge order the system remains conducting. 
At 88 K a magnetic phase transition occurs, associated 
with a structural distortion, and the onset of a long range 
magnetic order. Again, despite this transition, the sys- 
tem remains conducting up to 53 K where a small charge 
gap (14 meV) finally opens^ii^. The understanding of 
this complex phase diagram is still under process. Dif- 
ferent hypotheses have been advanced such as successive 
Fermi surface nestingil, polarons formationi^, one band 
versus three bands cross-over—. 



As already mentioned, the C0O2 compound remains 
conducting for all temperatures, opposite to the Mott in- 
sulator behavior that was expected due to its half-filled 
character. The low temperature dependence of the mag- 
netic susceptibility lead some authors to suggest that this 
compound is a weakly correlated metal^. However, as 
for the other compounds with x < 0.5, these data were 
reinterpreted^ in view of NMR measurements as due to 
strong antiferromagnetic coupling. 



In order to understand the low energy physics of these 
systems it is thus necessary to have accurate determina- 
tions of the pertinent degrees of freedom and of their cou- 
pling, as a function of the exact crystallographic struc- 
ture. Indeed, all authors agree on the importance of the 
sodium ordering and the induced cobalt local environ- 
ment distortions, on the low energy properties. The aim 
of this paper is to provide such informations using state 
of the art quantum chemical spectroscopy methods for 
the a; = 0.5 and x — compounds. The next section 
will shortly recall the methods. Section III will presents 
our results and finally the last section will sum up our 
conclusions. 



II. METHOD AND COMPUTATIONAL 
DETAILS 

A. Ab initio calculations 

The method used in this work (CAS+DDCpO) is a 
configurations interaction method, that is an exact di- 
agonalization method within a selected set of Slater 
determinants, on embedded crystal fragments. This 
method has been specifically designed to accurately treat 
strongly correlated systems, for which there is no single- 
determinant description. The main point is to treat ex- 
actly all correlation effects and exchange effects within a 
selected set of orbitals (here the 3d shell of the cobalt 
atoms) as well as the excitations responsible for the 
screening effects on the exchange, repulsion, hopping, etc. 
integrals. 

The CAS-I-DDCI method has proved very efficient to 
compute, within experimental accuracy, the local interac- 
tions (orbital energies, atomic excitations, exchange and 
transfer integrals, coulomb repulsion etc.) of a large fam- 
ily of strongly correlated systems such as high Tc cop- 
per oxides^'i, vanadium oxides^^, nickel and cuprate flu- 
orides^^, spin chains and ladders^**, etc. 

The clusters used in this work involve one and two 
cobalt atoms (CoOg and C02O10) and their oxygen first 
coordination shell (see figure [5]) . These fragments are 
embedded in a bath designed so that to reproduce on 
them the main effects of the rest of the crystal ; that is 
the Madelung potential and the exclusion effects of the 
electrons of the other atoms of the crystal on the clusters 
electrons. 

The electrostatic potential is reproduced by a set of 
point charges located at the atomic positions. The 
charges are renormalized next to the bath borders, us- 
ing the Gelle's method, in order to obtain an exponential 
convergence of the Madelung potential--. The conver- 
gence accuracy was set in the present work to the mili- 
electron-Volt. The nominal atomic charges used in this 
work are the formal charges, that is 4-3.5 and +4 for 
the cobalt atoms of the a; = 0.5 and x — compounds, 
—2 for the oxygen atoms and -t-1 for the sodium atoms. 
The exclusion effects are treated using total ions pseudo- 
potentialsSS (TIP) on the first shell of atomic sites sur- 
rounding the clusters. 




FIG. 2: a) CoOs and b) C02O10 fragments used in the present 
calculations. 

The calculations were done using the M0LCAS2i and 
CASDI^* suites of programs. The basis sets used can be 
found in reference |29| . The structural parameters were 
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taken for the Nao.5Co02 system from the 10 K data 
of reference [l^ and for the NaoCo02 system from ref- 
erence S 
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with energy 
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Extracting effective interactions from ab initio 
results 



The results of the ab initio calculations on an embed- 
ded fragment of the crystal provide us with energies (En) 
and eigenfunctions (4'„) for the local states from which 
the crystal ground and low lying excited states are built. 
Even though these eigenfunctions are expanded over mil- 
lions of determinants, they are dominated by only few of 
them, that arc usually considered as the implicit basis of 
the effective models (t — J or Hubbard in the present 
case). From the effective Hamiltonian theory (for in- 
stance such as derived in the quasi-degenerate pertur- 
bation theory^) one can extract the main concept nec- 
essary to get the effective Hamiltonian parameters from 
the ab initio results. Indeed, the effective Hamiltonian 
on the fragment should verify at the best the following 
set of equations 



^'2 

with energy 



\egiegieg2aigaig) 



2U + 8V- AJh 



(1) 



where En and '^n are the ab initio energies and eigen- 
states, and V is the projection operator on the local con- 
figurations supporting -ffcff- For instance, in the present 
case, these configurations can be the determinants sup- 
porting the local singlet and triplet state in the t — J or 
Hubbard model on the aig orbitals. The iJoff parame- 
ters are thus obtained from a mean-square fit minimiza- 
tion of equation [T] These mean-square fit equations can 
sometime be inversed and the parameters obtained from 
simple algebraic equations. This is for instance the case 
for the effective exchange J that is the singlet-triplet ex- 
citation energy of a C02O10 fragment. Similarly, the ef- 
fective transfer can be simply expressed (see equation [5] 
in section llH Cp from the same excitation energy and the 
knowledge of the wave-functions. 

One should however remember that the values of the 
effective integrals strongly depend on the nature of the 
effective model. Indeed the processes renormalizing the 
integrals in one model can be explicit in another and one 
should take this into account in order to avoid double 
counting. Let us illustrate our purpose on the t2g orbital 
energy splitting. In a Hubbard model based on the cobalt 
3d orbitals the CoOg fragment low lying states can be 



where U is the repulsion of two electrons on the same 
orbital, V and Jh the repulsion and exchange of two 
electrons on different 3d orbitals. It results that if Eg, 
El and E2 are the ab initio energies associated with these 
states, the ab initio spectrum should be associated with 
the orbital energies differences 
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and Eo 



Eo 



On the contrary if one extends the model in order to 
explicitely treat the differences — according to the na- 
ture of the 3d orbitals — in the two-electrons two-orbitals 
repulsion and exchange integrals, then the ab initio spec- 
trum does no more correspond to pure orbital energy 
differences but rather to 

Ei-Eq ^ Ea.g - Se^^ + 2(V„i,e^, ~ V^'^^ e',J 
~{Jh aig - Jh e'g^ e'^^ ) 

E2-E0 = eaig-e,,.^+2{V,,^,,^^V,,^,>J 

-{Jh aig e'g^ - Jh e'g^ e^^J 

In these equations the two electrons part is non zero as 
soon as the e'g orbitals are not pure t2g ones, but rather 
hg-Eg hybrids, as it is the case as soon as the trigonal 
distortion sets in (see reference [H and section IIIIBI for 
the present results). Whether the Racah's parameters 
are explicitely included in the model or the spherical ap- 
proximation used, one should thus use different Eq, — Eg' . 

9 gi 

effective orbital energy differences in order for the effec- 
tive Hamiltonian to fit the same ab initio spectrum. In 
the case of C0O2 layers it means different ordering for 
the aig and e'g orbitals (see reference IT^ . As usually 
done, we will use the spherical approximation in this work 
and renormalize the orbital energies by the effect of the 
Racah's parameters. 



III. RESULTS 



A. On site orbital energy splitting and charge 
order 



A strong controversy has been going on in the liter- 
ature about the cobalt 3d orbital splitting. Indeed, au- 
thors did not agree on the relative order of the t2g orbitals 
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under the trigonal distortion of the regular octahedron 
{T2g — >■ Aig (B Eg). The origin of the controversy was 
the fact that the crystalline field theorj^ as well as some 
LDA calculations'^^ found the aig orbital of lower energy 
than the e'g ones, while other density functional calcu- 
lations^, as well as quantum chemical calculations'^^ or 
ARPES experimental results^^ found them of reverse or- 
der. We recently showed^^ that the origin of the con- 
troversy was an incorrect treatment of the exchange and 
correlation integrals within the 3d shell. Indeed, only 
when these effects are taken into account with their whole 
complexity — that is when the orbital dependence of the 
two-electron two-orbital on-site repulsion and exchange 
are correctly treated — , the relative order between the 
aig and e'g orbitals issued from the t2g is correctly pre- 
dicted. Unfortunately this is not the case in the LDA ap- 
proximation (essentially due to the local exchange) and 
thus LDA predicts incorrect aig-e'^ ordering. This split- 
ting is however one of the crucial parameter for the de- 
termination of the nature of Fermi level states as recently 
shown by Marianetti et a^. Indeed, they showed (using 
DMFT calculations) that the existence of the e'g pock- 
ets in the Fermi surface of the NaxCo02 compounds are 
directly related to the this splitting. One can thus con- 
clude that the discrepancy between LDA and ARPES 
Fermi surfaces can be directly related to the improper 
treatment of the b Racah parameter in LDA. In fact, one 
can expect that an orbital dependant LDA-I-U, where the 
whole complexity of the Racah parameters would be in- 
cluded, should yield Fermi surfaces in agreement with the 
experimental ones. 

We thus computed the t2g orbitals splitting in the two 
compounds, for the different cobalt sites. These values 
can be extracted from the spectroscopy of the CoOe em- 
bedded fragments. Indeed, the excitation energy between 
the three cobalt states, where one hole is located on one 
of the t2g orbitals, yield the effective splittings between 
the aig and e'g orbitals. The relative order of the t2g 
orbitals is displayed on figure [3] One sees immediately 
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FIG. 3: Schematic representation of the effective cobalt t2g 
orbital energies in meV. Let us recall that the Co(l) cobalt of 
the a; = 0.5 system is located directly under or on top a Na^ 
ion. 

that the aig orbitals are always higher than the e'g ones 
as can be expected from reference 1 12l. In the a; = 0.5 sys- 
tem, both sites exhibit aig orbitals about 250-300 meV 
higher than the e'g ones. Similarly, in our dimer calcu- 



lations we see that the Co^+(l) ion is globally of lower 
energy than the Co''"'" (2) ion. Indeed, our calculations 
yield a global charge order of about 0.17 electrons in fa- 
vor of the Co(l) site. This result is in agreement with 
the common sense expectations that will favor a smaller 
valence for the cobalt with the largest number of Na"*" 
cations neighbors. Indeed, computing the electrostatic 
potential difference between the two cobalt sites one finds 
the Co(l) site about 315 meV lower than the Co (2) site. 
This result is however in contradiction with the Bond Va- 
lence Sum (BVS) model, since the later yields a charge 
ordering of 0.12 electron in favor of the Co(2) site^ how- 
ever it is in agreement with the LDA-I-U calculations^ 
(^ 0.16e in favor of Co(l)). In the pure C0O2 system, 
the aig-e'g orbital splitting is very large, much larger than 
in any of the other NaxCo02 systems, even larger than 
in the superconducting compound. This very large split- 
ting can be easily explained from (i) the large trigonal 
distortion (the largest observed in this family) and (ii) 
the very small Co-O distances. Indeed, it was shown in 
reference [12! that (i) the orbital splitting increases with 
increasing trigonal distortion and (ii) that the slope of the 
splitting as a function of the distortion increases linearly 
with decreasing Co-O distances. 



B. Orbital hybridization 

Let us now focus on the cobalt 3d orbital hybridiza- 
tion. These orbitals can hybridize in two ways : with the 
oxygen ligands 2p orbitals, but also within the cobalt 3d 
shell between the t2g and eg sets of the regular octahe- 
dron. Indeed, it was shown (see ref. [l^ for details) that 
not only this hybridization is symmetry allowed by the 
trigonal distortion of the octahedron, but also that it is 
responsible for the lowering of the e'g orbitals compared 
to the aig one. The t2g~eg hybridization angle is reported 
in table U and found to be non negligible, of the same or- 
der of magnitude as found in the other Co02-based sys- 
tems. The t2g-eg hybridization is larger for a: = than 
for X — 0.5, as can be expected from the larger trigonal 
distortion. 



Orbital 


Nao.5Co02 


C0O2 




Co(l) Co(2) 






11.59 11.54 


11.86 


692 


8.21 10.18 


11.86 



TABLE I: t2g-eg mixing angle (degrees) in the occupied e'g 
cobalt orbitals. 

As far as the cobalt ligand hybridization is concerned, 
we found it to be weak for the aig and e'g orbitals, with 
less than 10% weight on the oxygens. It is however quite 
large for the empty eg orbitals with 40 to 45 % weight on 
the neighboring oxygens. 
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C. Interatomic interactions 

Let us now focus on the interactions between two first 
neighbor cobalt atoms. Many authors assume that the 
low energy physics of this compound is determined by the 
aig orbitals only. The large aig-e'g orbital splitting found 
in our calculations seem to support this assumption. We 
thus computed the effective integrals both for a t — J 
model and a one-band Hubbard model. 

Analyzing the low temperature crystallographic struc- 
ture of Nao.5Co02 (10 K data of reference Hj), one sees 
that there are four types of independent Co-Co bonds 
(see figure SI), the Co(l)-Co(l) bond, the Co(2)-Co(2) 
bond and two kinds of Co(l)-Co(2) bonds, namely those 
where the cobalt-cobalt interactions are mediated by one 
0(1) and one 0(3) oxygen ligands and those where the 
cobalt-cobalt interactions are mediated by one 0(2) and 
one 0(3) oxygen ligands (see figure H]). The effective 




FIG. 4: (color online) Schematic representation of the C0O2 
layers in the Nao.5Co02 compound. There are two differ- 
ent cobalt sites represented in light and darker turkoise color, 
three different oxygen ligands and the four different Co- 
Co bonds represented in orange [Co(l)-Co(l)], red [Co(2)- 
Co(2)], dark blue [Co(l)-Co(2) bridged by 0(2) and 0(3) 
oxygens] and black [Co(l)-Co(2) bridged by 0(1) and 0(3) 
oxygens] . 

magnetic exchange integrals are 

J(ll) = -llmeV J(22) = -27meV 
j(12[i3]) = -i9meV J{12^^^^) = -19meV 

where J{ij) refers to the Co(i)-Co(j) bond and the [ab] 
superscript refers to the bridging oxygens. For the pure 
C0O2 system the magnetic coupling was computed to be 
much larger with 

J = -52meV 

In our notations, the negative integrals correspond to an- 
tiferromagnetic coupling. We thus find all exchange inte- 



grals, both in the a; = and x = 0.5 compounds, antifer- 
romagnetic, in agreement with the neutron scattering— 
and NMR-d0 data. One should notice that the exchange 
integrals are strongly modulated according to the local 
environment. Indeed, they differ by a factor larger than 
two for the different bond of the x = 0.5 system. It results 
that this system will be described by a strongly inhomo- 
geneous t — J model. In the x = compound the AFM 
exchange is very large, in agreement with NMR exper- 
imental results^. The difference in magnitude between 
the antiferromagnetic coupling of the x = 0, x = 0.5 and 
superconducting system'^ can easily be explained by the 
variations in the Co-0 distances. Indeed, these distances 
control for a large part the super-exchange term of the 
coupling that increases with decreasing distances. It re- 
sults that following De Vaulx et a? we found that AFM 
fluctuations increase with decreasing x. 

Let us now focus on the hopping integrals. One can 
get them either from the charge fluctuation in the singlet 
state of the Co^"*"- Co"^"*" fragments, or from the spec- 
troscopy of the Co^"*"- Co*"*" fragments. A simple analy- 
sis shows that the two hopping integrals are somewhat 
different. In the first case, there is only one spectator 
electron on the bond, while in the second case, there are 
two of them (see figure [5]). It results in a different eval- 
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FIG. 5: Schematic representation of the effective hopping in- 
tegrals as a function of the spectator electrons on the bond. 
I and r labels refers to the arbitrary "left" and "right" atoms 
in the fragment. 



nation of the hopping integrals according to the number 
of spectator electrons. 

t(le) = {di\h\dr) + {didi\l/ri2\drdi) 

t{2e) = {di\h\dr) + {didi\l/ri2\drdi) + {didr\l/ri2\drdr) 

where the I and r labels refers to the arbitrary "left" and 
"right" atoms in the fragment, and h is the single electron 
Hamiltonian part. 

t{le) can be extracted from the Co''+-Co''+ fragment. 
Indeed, mapping a single band Hubbard model on the 
computed low-energy singlet and triplet states wave func- 
tions and energies, one gets after a little algebra 



t{le) = J 



1 



2 cos (p tan a 



U-V = J l- 



1 



cos 2(p tan a 



(2) 
(3) 



assuming the following form for the computed singlet and 
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triplet wave functions 

\didr) + \drdi) 



cosa- 



sma ^cos(^- 
small terms 
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mediated by the oxygen ligands) between the 3di orbital 
of one cobalt and idj of the other, tpij are effective intra- 
atomic transfer integrals between the di and dj orbitals 
of the same atom. One may be surprised to find such 
terms, since the direct integrals are zero due to the 
symmetry, however the coupling with the bridging oxy- 
gen 2p orbitals yield in second order perturbation theory 
an effective intra-atomic transfer term of 



\didr 



\drdi) 



V2 



small terms 



J = Esg — Etp- U — V is the effective on-site repulsion of 
a simple Hubbard model. Let us remind that it accounts 
for both the average on-site repulsion between the two 
sites, U = {Ui + Ur) /2, and the nearest neighbor effec- 
tive repulsion between the cobalt atoms, V — Vw- The 
(electron-)hopping integrals between the aig orbitals thus 
come as 



iiO(ll)=84meV tll{22) 



•-oov 

115 meV ioo(12[^^l) 



139 meV 
: lllmeV 



for the X — 0.5 compound and as 



I'qI = 212 meV 

for the X = one. 

Once again, one sees that the effective transfer inte- 
grals are strongly dependant on the local environment of 
the cobalts and that the Nao.5Co02 system is more inho- 
mogeneous than first thought. In the pure compound, the 
effective transfer between the aig orbitals of NN atoms 
is larger than both in the x = 0.5 and the superconduct- 
ing system. The relative variations of the Co-0 and large 
Co-Co distances between these three systems support the 
idea that the effective transfer integrals are dominated by 
the ligand-mediated term rather than the direct hopping 
term, despite the short metal-metal distances in these 
edge-sharing compounds. 

The transfer integrals with two spectator electrons on 
the bond, i^*^, can be extracted from the low-energy states 
of the Co'^"'"- Co^"*" fragments. There are six low energy 
states according to the localization of the hole on the aig 
and e'g orbitals. The associated effective Hamiltonian can 
be written as 
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(4) 



where are the atomic effective orbital energies (see fig- 
ure [3]), tfj are the effective transfer integrals (direct plus 



p 



where i and j refers to the id orbitals of the same cobalt 
atom, the sum over p runs over all the ligand bridging 
orbitals, tip is the cobalt 3di-ligand hopping integrals 
and Ap is the excitation energy toward the ligand-to- 
metal charge transfer configuration. For a more detailed 
description of the underlying mechanism, one can refer 
to reference [sli. Summing the tpij contributions coming 
from the six oxygens around a cobalt atom, one can show 
that they exactly cancel out in a symmetric system. In 
the present system, where the atomic S'g symmetry is not 
exactly respected, these terms will certainly not exactly 
cancel out, however, one can expect that their sum will 
remain very weak. 

Table |TT] displays the different effective integrals in- 
volved in matrix for the x = 0.5 system. One sees 



Bond 



Co(l)-Co(l) 
Co(2)-Co(2) 



Inter-atomic transfers (meV) 



,2e 
^00 



,2e 
'22 



,2e 



,2g 
I20 



225 -261 -4 -22 -18 -27 
281 -325 -11 ± 10 -18 ± 10 



TABLE II: Effective (electron-)hopping integrals between aig 
and e'g orbitals of two neighboring cobalt atoms (in meV). 

immediately that the dominant transfer integrals are the 
inter-atomic hopping between two Oi^ orbitals and two 
e'gi orbitals, in agreement with what we found on the 
superconducting compound. All other integrals are rela- 
tively weak (< 30meV) and can probably be omitted in 
a simple model. Comparing these values with what we 
found on the superconductor system, one sees that they 
are of the same order of magnitude, except for the inter- 
atomic hopping between the aig orbital of one cobalt 
and the e^j orbital of the other. Indeed, the larger value 
found in the superconducting systems (-53 meV) is now 
spread over the three ijg, t^Q and hopping terms, due 
to local symmetry breaking. 



D. On-site repulsion 

Let us now examine what we find for on-site the re- 
pulsion U from our calculations. Indeed, as detailed in 
equation [3l the value oiU —V can be extracted from the 
singlet and triplet energies and wave functions on the 
00"*+- Co^+ fragments. Results for a single band model 
are presented in table IIIIl One sees immediately that 
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Compound 




ITT ■\/\ U 

(U - V)/t 


C0O2 


3.71 


15.6 


[Nao.35Co02][H20]l.3 


3.65 


19.4 


Nao.5Co02 






[/(!)- 


2.61 


31.1 


(7(2) - y(22) 


2.86 


20.6 


f7- V(12["l) 


2.76 


24.0 


17-1/(12^^1) 


2.56 


23.1 



TABLE III: Effective on-site Hubbard repulsion (eV) and cor- 
relation strength for different NaxCo02 compounds. U = 
[U{l) + U{2)]/2 



the Hubbard on-site repulsion decreases with increasing 
sodium content, that is with increasing number of elec- 
trons in the C0O2 layers. Once again correlating this fact 
with the Co-0 distances one sees that the larger the Co- 
O distances, the weaker the cobalt on-site repulsion. This 
can be understood by the fact that a larger volume of the 
coordination shell around the metal atom will allow its 
magnetic orbitals to spread over some more and this re- 
duce the intra-orbital repulsion. If one now looks at the 
correlation strength (U/t) rather than the on-site repul- 
sion value, one sees that the tendency is reverse. Indeed, 
the larger the sodium content, the larger the correlation 
strength. It is easy to see that the opposite variation of 
the correlation strength and the one-site repulsion, as a 
function of x, is again due to the Co-0 distances. Indeed, 
if an increase of the metal coordination sphere induce a 
decrease of the effective U, it induces a much stronger 
reduction of the effective hopping integral t, thus result- 
ing in weaker electronic correlation. This results fit quite 
well with what was originally supposed (from experimen- 
tal observations) in these systems, that is a decrease of 
the correlation strength with x. However it seems hard 
to say that these systems are weakly correlated, since the 
smallest U/t ratio is 15. 



IV. CONCLUSION 

We determined the effective on-site and nearest- 
neighbor coupling parameters for the the Nao.5Co02 and 
C0O2 compounds within Hubbard and t — J models. The 
effective integrals and orbital energies were computed us- 
ing ab-initio quantum chemical methods treating exactly 
the strong correlation effects within the cobalt 3d shell 
and the screening effects up to the double-excitations. 

The Hubbard on-site repulsion and the correlation 
strength U/t were found to vary in opposite ways with 
the sodium content. Indeed, when doping the C0O2 
layer, the t value decreases more rapidly than the U value 
and the resulting correlation strength increases. For the 
X = limit, however, the system is still strongly cor- 
related with U/t ^ 15. This result is in disagreement 



with Onoda's conclusion that the C0O2 compound is a 
weakly correlated metal^ii^. On the contrary, Julien^ 
concludes the system is a strongly correlated metal close 
to the Mott phase transition. Our results clearly do 
agree with the latter hypothesis, indeed, it was shown 
using DMFT calculation that the Hubbard model reach 
the Mott transition in a triangular lattice at a critical 
correlation strength of Uc/tc = 15^. Our evaluation of 
the correlation strength : 15.6, clearly situates the C0O2 
compound very close to the Mott transition. Similarly, 
our results are in agreement with NMR measurements^ii^ 
that sees an increase of the antiferromagnetic coupling 
with decreasing sodium content. However, contrarily to 
Kawasaki et we do not find the pure C0O2 compound 
behaving on a specific way. Indeed, we found it to dis- 
play the largest antiferromagnetic coupling among all the 
systems we studied. 

We determined the ligand field splitting for the dif- 
ferent crystallographically independent cobalts of the 
Nao.5Co02 and C0O2 compounds. As for the super- 
conducting compound, we find that the aig orbitals are 
of higher energy than the e'g ones, yielding a e'g aig^ 
cobalt atomic configuration. Comparing the t2g split- 
ting in these two compounds with the splitting in the 
superconducting compound, we see that it follows the 
behavior predicted in our previous work [l^ both as a 
function of the amplitude of the trigonal distortion and 
as the size of the cobalt first coordination shell. Indeed, 
in the C0O2 compound, the average distance between 
the oxygen and cobalt planes is 0.85A while it is 0.88A 
in the superconducting system and ranges between 0.95 
and 0.98A in the x — 0.5 compound. As this distance 
increases with increasing electron doping of the C0O2 
layer, the splitting between the e'g and aig orbitals de- 
creases, being maximum for the pure system. Comparing 
this orbital order with the LDA or LDA-I-U calculations 
one sees that it is reverse order. Indeed, tight-binding 
fitting of the LDA+U calculations yield a higher hole en- 
ergy by about O.OleV when located in the aig orbital^. 
Why do LDA calculations yield a reverse order between 
the aig and e'g atomic orbitals than wave-function clus- 
ter calculations? The answer is in the treatment of the 
Coulomb and exchange terms within the 3d shell. In- 
deed, in LDA calculations the correlation and exchange 
contributions are treated using a local approximation. It 
means that the small differences in energy that should 
occur in these term according to the nature of the oc- 
cupied 3d orbitals (the b Racah's parameter) is ignored. 
This is even worse with LDA-f U calculations since the 
Hartree term is also replaced and a unique U (or rather 
U, V and Jh) is used while the Coulomb and exchange 
integrals between two different 3d orbitals depend on the 
nature of the considered orbitals. These integrals differ- 
ences are quite small and, in most cases, the averaging 
is valid enough. In the present case, however, it can 
be shown that the 2-electron contribution to the orbital 
energy that is responsible for the aig-e'g inversion in ref- 
erence is a direct function of the Racah's b integral, 
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and of the Coulomb and exchange integrals differences. 
The effects of these imperfect exchange and correlation 
treatment within the 3d shell, and thus of the resulting 
orbital order between the aig and orbitals is to pop 
up the band issued from the later in comparison with 
the bands issued from the former. Indeed, it was shown 
by Marianetti et (using DMFT calculations) that 
when the aig orbital is supposed of higher energy (lower 
hole energy) than the orbital, the e'g pockets in the 
Fermi surface disappear. More recently Bourgeois et 
showed that using orbital splitting of the same sign and 
order of magnitude as in the present work, DMFT Fermi 
surfaces and band structure nicely fit the ARPES data 
(they studied different systems ranging from x = 0.3 to 
X = 0.7). In view of our results and the above analysis we 
think that the e'g Fermi pockets found in LDA calcula- 
tions and not observed in ARPES experiments are due to 
an intrinsic feature of the LDA exchange and correlation 
functional and thus that these pockets are true artefacts 
of the method. 



Finally, one should remember the strong variation of 
the effective exchange and hopping integrals as a function 
of the local distortions. This fact can be of importance in 
the understanding of the changes in magnetic behavior, 
as a function of x, in the NaxCo02 family. It would thus 
be of interest to see if the local distortions observed for 
larger x values will be able to reverse the sign of the 
magnetic exchange as suggested by neutron diffraction 
experiments. 
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